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Abstract — One of the popular dynamics on complex networks 
is the epidemic spreading. An epidemic model describes how 
infections spread throughout a network. Among the compartmen- 
tal models used to describe epidemics, the Susceptible-Infected- 
Susceptible (SIS) model has been widely used. In the SIS model, 
each node can be susceptible, become infected with a given 
infection rate, and become again susceptible with a given curing 
rate. In this paper, we add a new compartment to the classic 
SIS model to account for human response to epidemic spread. 
Each individual can be infected, susceptible, or alert. Susceptible 
individuals can become alert with an alerting rate if infected 
individuals exist in their neighborhood. An individual in the alert 
state is less probable to become infected than an individual in 
the susceptible state; due to a newly adopted cautious behavior. 
The problem is formulated as a continuous-time Markov process 
on a general static graph and then modeled into a set of 
ordinary differential equations using mean field approximation 
method and the corresponding Kolmogorov forward equations. 
The model is then studied using results from algebraic graph 
theory and center manifold theorem. We analytically show that 
our model exhibits two distinct thresholds in the dynamics of 
epidemic spread. Below the first threshold, infection dies out 
exponentially. Beyond the second threshold, infection persists 
in the steady state. Between the two thresholds, the infection 
spreads at the first stage but then dies out asymptotically as the 
result of increased alertness in the network. Finally, simulations 
are provided to support our findings. Our results suggest that 
alertness can be considered as a strategy of controlling the 
epidemics which propose multiple potential areas of applications, 
from infectious diseases mitigations to malware impact reduction. 



I. Introduction 

Modeling human reactions to the spread of infectious dis- 
ease is an important topic in current epidemiology IT], 121, 
and has recently attracted a substantial attention lf3l- lfToll . 
However, few papers are available in the literature which 
consider the human response to the epidemic in a systematic 
framework and the contributions to the problem are still in 
an early stage. The challenges in this topic concern not only 
how to model human reactions to the presence of epidemics, 
but also how these reactions affect the spread of the disease 
itself. In a general view, human response to an epidemic 
spread can be categorized in the following three types: 1) 
Change in the system state. For example, in a vaccination 
scenario individuals go directly from susceptible state to 
recovered without going through infected state. 2) Change 
in system parameters. For example, as in LIU , individuals 
might choose to use masks. Those who use masks have a 
smaller infection rate parameter, 3) Change in the contact 
topology. For example, due to the perception of a serious 
danger, individuals reduce their contacts with other people who 
can potentially be infectious ^. 



Early results on epidemic modeling dates back to lfT2l . In 
ifTJi an epidemic model on a homogenous network was stud- 
ied. Later on, results for heterogeneous networks were reported 
in lfT4l . Pastor-Satorras et. al. JTSI studied epidemic spreading 
in scale free networks, showing that in these networks the 
epidemic threshold vanishes with consequent concerns for the 
robustness of many real complex systems. Wang et. al. llT6l 
provided the first result for a non-synthetic contact topology, 
and studied the epidemic spread dynamic on a general static 
graph. Through a local analysis of a mean-field discrete model, 
it was shown that the epidemic threshold is directly related 
to the inverse of the spectral radius of the adjacency matrix 
of the contact graph. More detailed proof was provided in 
fVi\. Ganash et. al. |18| proved the same result without 
any mean-field approximations. A continuous-time epidemic 
model was studied by Van Mieghem et. al. ||T9l . where a set 
of ordinary differential equations was extracted through mean- 
field approximation of a continuous time Markov process. 
The relation between the epidemic threshold and the spectral 
radius was rigorously proved and further insights about the 
steady state infection probabilities were analytically derived. 
Preciado and Jadbabaie ll20l studied the epidemic spread on 
geometric random networks and then in f2T1, they investigated 
the epidemic threshold on a general contact graph with respect 
to the network structural information. 

A good review on existing results in the literature where 
the human behavior is taken into account for epidemic mod- 
eling can be found in |2). Poletti et. al. 1221 developed a 
population-based model where susceptible individuals could 
choose between two behaviors in response to presence of 
infection. Funk et. al. JH showed that awareness of individuals 
about the presence of a disease can help reducing the size 
of the epidemic outbreak. In their paper, awareness and 
disease have interconnected dynamics. Theodorakopoulos et 
al. f3\ formulated the problem so that individuals could make 
decision based on the perception of the epidemic size. Most 
of the existing results are suitable for a society of well-mixed 
individuals, since the contact graph is usually considered to 
be homogeneous (i.e. all nodes have the same degree). To 
the authors' knowledge, the study of the human response in a 
realistic network of individuals with a general contact graph 
has not been reported so far 

In this paper, we model the human response to epidemic 
in the following way. A new compartment is considered in 
addition to susceptible and infected states. A susceptible indi- 
vidual becomes alert with some probability rate if surrounded 
by infected individuals. An alert node gets infected with a 
lower rate compared to a susceptible node does with the 
same number of infected neighbors. The contribution of this 



paper is two-fold. 1) Unlike most of the previous results, 
no homogeneity assumption is made on the contact network 
and the human-disease interaction in this paper is modeled 
on a general contact graph. 2) We show through analytical 
approaches that two distinct thresholds exist. The two are 
explicitly computed. To the authors' knowledge the existence 
of two distinct thresholds is reported for the first time in this 
paper, providing a fundamental progress on previous results. 
Additionally, this result has the potential to be applied to 
mitigate epidemics in several different complex systems, from 
human and animal infectious diseases, to malware propagation 
in computer and sensor networks. 

The rest of the paper is organized as follows. In Section HH 
some backgrounds on graph theory, center manifold method, 
and the N-Intertwined SIS model (developed in [19] ) are 
recalled. Section |III] is devoted to the problem formulation and 
model derivations. Stability analysis results of the model are 
provided in Section |IV] Finally, results are examined through 
numerical simulations in Section IVl 

II. Preliminarily and Background 

A. Graph Theory 

Graph theory (see |23|) is widely used for representing the 
contact topology in an epidemic network. Let Q — {V,£} 
represent a directed graph, and V = {1, ..., A^} denote the set 
of vertices. Every individual is represented by a vertex. The 
set of edges is denoted as £ C V x V. An edge is an ordered 
pair {i,j) G £ if individual j can be directly infected from 
individual i. In this paper, we assume that there is no self 
loop in the graph, that is, (i, i) ^ £. Mi = {j G V \{j, i) G £} 
denotes the neighborhood set of vertex i. Graph Q is said to 
be undirected if for any edge {i,j) G £, edge {j,i) E £. A 
path is referred by the sequence of its vertices. A path P of 
length k between vq, Vk is the sequence {vo,...,Vk} where 
{vi-i,Vi) G £ for i = l,...,fc. Directed graph Q is strongly 
connected if any two vertices are linked with a path in Q. 
A = [ttij] G R^^^ denotes the adjacency matrix of Q, where 
Uij = 1 if and only if {i,j) G £ else Uij — 0. The largest 
magnitude of the eigenvalues of adjacency matrix A is called 
spectral radius of A and is denoted by p{A). 



B. Center Manifold Theory 

Linearization is a useful technique for local stability analysis 
of nonlinear systems. However, in the cases where lineariza- 
tion results in a linear system with some negative real part and 
some zero real part eigenvalues, the linearization method fails. 
In these cases, the local stability analysis can be performed by 
analyzing a nonlinear system of the order exactly equal to the 
number of eigenvalues with zero real parts. This method is 
known as center manifold method. In this section, we have a 
quick review on center manifold theory. More details can be 
found in 1 24 1 and ll25l . 

For Zs G M"^ and Zc G K"% consider the following system 



Zs — -AgZg ~\- gs{ZcT Zg) 



(1) 

(2) 



where the eigenvalues of Ag G R"==<"» and Ac G ]R"<^=<"<^ 
have negative and zero real parts, respectively. The functions 
gc and gs are twice continuously differentiable and satisfy the 
conditions 



ffi(0, 0) = 0, V.g,(0, 0) = 0, i G {s, c}, 



(3) 



where is a vector or matrix of zeros with appropriate 
dimensions. There exists a function h : M"= -^ M"= satisfying 



h{Q) = 0, V/i(0) = 0, 



(4) 



that Zs = h{zc) is an invariant manifold (see ll25l for the 
definition) for ([T]i and dU near the origin. The dynamic system 
dU and (|2|i can be studied through the reduced system 



Zc = AcZc+gc{Zc,HZc))- 



(5) 



The invariant manifold Zs = h{zc) is a center manifold for 
the system ([T]i and (|2]l, i.e., every trajectory of ^ and (|2|l 
with the initial condition Zc(0) = Zc{0) and ^^(0) — h{zc{0)) 
satisfies Zc{t) = Zc{t) and Zs{t) ~ h{zc{t)). In addition, 
small deviation from the center manifold is exponentially 
attracted, i.e., if ||zs(0) — /i(ic(0))|| is small enough, then 
||zs(i) — h{zc{t))\\ will go to zero exponentially. 

C. N-Intertwined SIS Model for Epidemic Spread 

We have built our modeling based on a newly proposed 
continuous-time model for epidemic spread on a graph. Van 
Mieghem et. al. |19| derived a set of ordinary differential 
equations, called the N-intertwined model, which represents 
the time evolution of the probability of infection for each 
individual. The only approximation for the N-intertwined 
model corresponds to the application of the mean-field theory. 

Consider a network of N individuals. Denote the infection 
probability of the i-th individual by pi G [0, 1]. Assume that 
the disease is characterized by infection rate /3q G K+ and cure 
rate 5 G M^. Furthermore, assume that the contact topology 
is represented by a static graph. The N-intertwined model 
proposed in |19| is 



where a^ == 1 if individual j is a neighbor of individual i, 
otherwise a^ = 0. 

Proposition 1: Consider the N-intertwined model (|6]l. Ini- 
tial infection will die out exponentially if the infection strength 

T ^ -f- satisfies 



.aA,< 



1 



(7) 



6 - p{Ay 

where p{A) is the spectral radius of the adjacency matrix A 
of the contact graph. 

Remark 1: The value r^, = —^ is usually referred to as 
the epidemic threshold. For any infection strength r > r^, 
infection will persist in the steady state. The following result 
discusses the steady state values for infection probabilities. 

Proposition 2: If the infection strength is above the epi- 
demic threshold, the steady state values of the infection 



probabilities, denoted by pf for the i-th individual, is the 
non-trivial solution of the following set of equations 



continuous-time Markov process: 



A) 
s 






PI 



te{i,...,N}. 



(8) 



III. Model Development 

In this paper, we add a new compartment to the classic SIS 
model for epidemic spread modeling to propose a Susceptible- 
Alert-Infected-Susceptible (SAIS) model. The contact topol- 
ogy in this formulation is considered as a general static graph. 
Each node of the graph represents an individual and a link 
between two nodes determines the contact between the two in- 
dividuals. Each node is allowed to be in one of the three states 
"5'.- susceptible", "I: infected", and "A: alert". A susceptible 
individual becomes infected by the infection rate /3q times 
the number of its infected neighbors. An infected individual 
recovers back to the susceptible state by the curing rate 6. An 
individual can observe the states of its neighbors. A susceptible 
individual might go to the alert state if surrounded by infected 
individuals. Specifically, a susceptible node becomes alert 
with the alerting rate k G M+ times the number of infected 
neighbors. An alert individual can get infected in a process 
similar to a susceptible individual but with a reduced infection 
rate < /3^ < /3q. We assume that transition from an alert 
individual to a susceptible state is much slower than other 
transitions. Hence, in our modeling setup, an alert individual 
never goes directly to the susceptible state. The compartmental 
transitions of a node with one single infected neighbor are 
depicted in Fig. [T] 




Fig. L The infected population fraction in Example, (a) SIS model, (b) SAIS 
model with, (c) SAIS model with . The Monte-Carlo simulation results are 
shown in blue. 



The epidemic spread dynamic is modeled as a continuous- 
time Markov process. For each node i G {1,...,A^}, define 
a random variable Xi : {S,I,A} -^ {0,1,2}. Denote X* 



Pr[Xr^* - l\X! - 0] = /3„At J2 l{x-i} + o{At), 

Pr[X*+^* = Q\Xl = 1] = (5Ai + o(At), 

Pr[X*+^* = 2\Xl = 0] = kM J2 1{x|=i} + o{At), 

Pr[X*+^* - l\Xl - 2] = P,M Y, l^x|=i} + o(At), (9) 

for i e {1, ..., N}. In Q, Pr[-] denotes probability, Ai > is 
a time step, and l{x} is one if X is true and zero otherwise. 
A function /(At) is said to be o(At) if limAt-^o ^^ = 0. 
A common approach for studying a continuous-time Markov 
process is to derive the corresponding Kolmogorov forward 
(backward) differential equations (see I26J and I27J ). As can 
be seen from the above equations, the conditional transition 
probabilities of a node are expressed in terms of the ac- 
tual state of its neighboring nodes. Therefore, each state of 
the Kolmogorov differential equations corresponding to the 
Markov process ^ will be the probability of being in a 
specific configuration. In this case, we will end up with a 
set of first order ordinary differential equations of the order 
3^. Hence, the analysis will become dramatically complicated 
as the network size grows. In addition, it is more desirable 
to study the probability that each individual is susceptible, 
infected, or alert. Using a proper mean-field approximation, 
it is possible to express the transition probabilities in terms 
of infection probabilities of the neighbors. Specifically, the 
term l{x'=i} is replaced with Pr[X* = 1] in dU. Hence, the 
following new stochastic process is obtained: 



Pr[X*+'^* = 1\XJ = 0] = /3oAt J2 Pi'l^j = 1] + o(^0, 

Pr[X*+^* = 0\Xl = 1] = SAt + o(At), 

Pr[X*+^* = 2|X* = 0] = nAt J^ ^^i^j = 1] + "(At), 

l\Xl = 2] = /3„At ^ Pr[XJ = 1] + o{At). 

(10) 



Pr[X 



t+At 



Define a new state Xi = [si,Pi,qi\^ , where Si, pi, and 
qi denote the probabilities of individual i to be susceptible, 
infected, and alert, respectively. The Kolmogorov forward 
differential equations of the stochastic process ( fTOb can now 
be found as 



QJx^, ie{l,...,A^}, 



(11) 



where 



e. 



~5 





6 


Pa E aijPj 
Pa E 0■^JPj 


-Pa E a.^JPj 




-iPo + '^) E o-ij?] 

(12) 



a measure of the random variable at time t for node i. is the infinitesimal ti-ansition matrix. One property of the 
The epidemic spread dynamics is modeled as the following dynamic system ( fTTI ) is that Si+pi + qi is a preserved quantity. 



Hence, the states Si, pi, and qi are not independent. Omitting Sj 
in ( fTTT i. the following set of differential equations is obtained: 



Pi = /3o(i -Pi- qi) XI "'J^J + ^a9» 



E 



ijjPj - 5p.i, 



k(1 -pi -gi) 



/ . "ij/ 



(13) 

foriG {1,...,A^}. 

Remark 2: As can be seen using a mean-field approxima- 
tion, the dimension of the differential equations is reduced 
from 3^ to 2N . However, some information is definitely lost 
and there is some error For example, the Markov process 
(|9]l exhibits an absorbing state. However, no absorbing state 
can be observed based on the equations ( fT3T l and (fl4] i. In 
addition, as is discussed in |fT9l . the solution from the mean 
field approximation is an upper-bound for the actual model. 

IV. Behavioral Study of SAIS Epidemic Spread 
Model 

In this section, the dynamic system ( fTsT l and (fl4] i derived 
in the previous section is analyzed. It is shown that alertness 
decreases the size of infection. In addition, in an SAIS epi- 
demic model, the response of the system can be categorized in 
three separate regions. These three regions are identified with 
two distinct thresholds t]. and r^. Below the first threshold, the 
epidemic dies out exponentially. Beyond the second threshold, 
the epidemic persists in the steady state. Between t]. and 
r^, the epidemic spreads at the first stage but then dies 
out asymptotically as the result of increased alertness in the 
network. 

A. Comparison between SAIS and SIS 

In this section, the SAIS model and the SIS model are com- 
pared in the sense of infection probabilities of the individuals. 
Specifically, we are interested to compare Pi{t), the response 
of ( fTSl l and (fl4] i. with infection probability p^(t) in the N- 
intertwined SIS model, which is the solution of the system 



P'l = /3o(l -P'<) y^ a,jp'. - 5p[. 






(15) 



It is shown that alertness decreases the probability of infection 
for each individual. This result is stated as the following 
theorem. 

Theorem 1: Starting with the same initial conditions 
Pi(io) = p'iito), i = {l:---:^}i the infection probabilities 
of individuals in SIS model ( fTsl l always dominate those of the 
SAIS model O and (0, i.e., 

P^{t)<p',{t),i^{l,...,N] Vie[io,H- (16) 
Proof: Rewrite the equations (fTSl l as 

Pi =^o(l"-PiO E "■vPj-{I^O-l^a)'it X a^jPJ-5p^. (17) 

Starting with the same initial conditions Pi(io) — p'iih), it is 
concluded that 

p,{to)^p',{to)^P;{to)<P-{to), 



since (3^ < /3q by definition and therefore (/3q — 
Pa)Qi{to) Z^jgM aijPjih) is a non-negative term. According 
to dlSl l. there exists tf > ig so that 

P^it)<p'M,^^{'^,■■■,N} yte[to,tf]. (19) 

The theorem is proved if we show that inequality (fT9] l holds 
for every tf G (io,oo). Assume that there exists ti > to, so 
that iT% holds for tf = ti but it is not true for any tf > ti. 
Obviously, at t = ii, 

3i e {1, ..., N} so that p,{ti) = p'^{ti) and p,{ti) > p'^iti). 

(20) 
In the subsequent arguments, it is shown that no such ti exists. 
From ( [TtI i. Pi{ti) is found to satisfy 

Pt{ti) = Pa{'^-Pi{ti))^ a',jpj{ti) 

< /3o(l --PK^i)) E "yPj(^i) ~ ^Piiti) 

= /3o(l -p^ii)) E ^^jPAti) - ^pU^i), (21) 

according to (l20l l and the fact that (/Jo^/^a)?*!^!) S '^ijPji^i) 
is a non-negative term. Based on (fT9] l. Vj G {!,..., A^} we 
have Pj{ti) < Pj(ii). Therefore, the inequality (ISTT i is further 
simplified as 

Mh) < M^-pKh)) E a^Jp',itl) - 6p',iti) ^p'Si)- 

(22) 
Having Pi{ti) < Pi{ti) contradicts with ( |20] |. Therefore, no 
such ii exists so that (|20] | is true. As a result the inequality 
(fT9] l holds for every tf G (io,oo). This completes the proof. 



B. Exponential Epidemic Die-Out 

Theorem 2: Consider the SAIS epidemic spread model ( fT3b 
and (fT4] i. Assume that the infection strength satisfies 



Po ^ 1 
5 P(A)- 



(23) 



(18) 



Then, initial infections will die out exponentially. 

Proof: The solution of pi (t) was proved in Theorem [1] 
to be upper-bounded by p^{t). As described in Section Hl-CI 
and based on Proposition [T] the N-intertwined model ( fTSl l is 
exponentially stable if ( l23T l is satisfied. As a consequence, Pi(i) 
in ( [Tt] ) is also exponentially stable if ( l23T l is satisfied. ■ 

Remark 3: In the proof of Theorem [T] g^ is considered as 
a non-negative time-varying term. Under the conditions of 
Theorem |2l qi (i) will regulate at some value depending on 
the initial conditions. 

Remark 4: Note that adding the alert compartment does not 
contribute to the epidemic threshold for exponential die out. 
This result is already concluded in |8| for a homogeneous 
network (i.e. all nodes have the same degree). 



C. Asymptotically Epidemic Die-Out 
According to (fl4] l. 



with 



ffi,,(p,r) ^ _{/3„ + ^2±^}p,^ 



<li = 



1 + 4- 






, *G{1,...,7V}, 



(24) 



-(/5o-/5a>^ 



is an equilibrium for ( fT4l i. To facilitate the subsequent analysis, 
define a new state r, as 



52,j(p,r) 






(35) 
(36) 



A e 1 - P^ 

! + £» 



(25) 



The derivatives pi and r^ in the new coordinate can be found 

by substituting qt = rt^ ^ 2^ from (|25] | in (fT3] l and 

(IT4] | as 



ft = /3o(l -ft - {''» + 



1 



Pj 



■})E 



11^ 11^ 



lijPj 



If we linearize the system (|30] | and dSTl i. the resulting system 
has N zero eigenvalues. Therefore, linearization technique 
fails to investigate the stability properties of (l30l l and (ISTT l. In 
the following arguments, we show that center manifold theory 
can be employed to study the stability of (l30l ) and JsTl i. 

The eigenvalues of matrix {P^gA — 61) are /3e„Ai — 6,i ^ 
{1, ...A^}, where A^'s are the eigenvalues of the adjacency 
matrix A. Therefore, assuming that 



+Pa{n 



Pa 



1 



1 



-^} ^ a,yPj -(5p,; 



/3 



eg 



< 



P(^)' 



(37) 



0-— 37 -1- A^a— — 97} XI "»J'PJ' 



-{/3o + 7^}ft E «^^-P^- 



(26) 



and 



the matrix {I3^„A — SI) is Hurwitz (i.e., a matrix that all of 
its eigenvalues have negative real parts). In addition, the two 
nonlinear functions Gi and G2 defined in ( l33T l and ( |34] | satisfy 

Gj(0,0) = 0, VGj(0,0) = 0, (38) 

for j e {1,2}. Hence, the center manifold theory reviewed in 
Section lTl-Bl may apply. The center manifold theorem suggests 
that there exists a function H{-) : M^ -^ M.^ where the 
dynamics (l30l l and ( |3TI ) can be determined by 



k(1 -p,- {r, + — - - — '—}) J2 "'^iPj 

1 + lf 1 + lf ]'-'' 

-Pain + —-j^ - 7-^} E °'J? 



f = G2(i/(f),f) 



(39) 






Differential equation ( |39] ) can be written in terms of its 
entries as 



K j<^Mi 



J6M 



JPJ- 



To facilitate the subsequent analysis, define 



(27) 



(28) 
(29) 









aij7ij(f), 



(40) 



According to ( |26] l and d27] ) and the definitions ( |28T l and ( |29] l, 
the followings are true 



for i £ {!,..., A^}, where /li(-) is the i-th component of 

Remark 5: Usually, it is not feasible to find hi(-) explicitly. 
In the subsequent analysis, instead of explicit calculations, we 
make use of the following property of hi{-): Since the prob- 
ability Pi is non-negative, each function hi{-) is necessarily 
non-negative. 

Lemma 1: The trajectories of ( l40b will asymptotically con- 
verge to the set defined by 



where 



and 



p = {P,gA-6I)p + Gi{p,r), 
r = Or + G2(p,r), 



Pen — PO , R ~^ Pa' 



1+^ "'^1+^^ 



Gi(-) ^ [5i,i(-),-,5i.m(-)]^, 
G2(-) ^ [92A-),...,92M-)f, 



(30) 
(31) 



(32) 



(33) 
(34) 



n = {re 



dW 



fi 2^ aijhjir) ==0}. 



(41) 



Proof: Define a continuously differentiable function V as 



V ^ if^f . 
2 



(42) 



Taking the derivative of V with respect to time, we have 



V 



N 

E ^»^'' 



-k(1 



i=l 



i=l 



f.2 sr 






ay/ij(f) . (43) 



It can be seen that the time derivative V is negative semi- 
definite according to Remark |5] According to the LaSalle's 
invariance theorem (see |25|) the trajectories of ( |40] i will 
asymptotically converge to the set V = 0, i.e., 



f] ^ {f e 



dW| 






aijhj{r) = 0}. 



(44) 



Theorem 3: Consider the SAIS epidemic model (fTsT i and 
( fT4l i. Assume that the infection strength satisfies ( [37] ) where 
(3 is defined in ( [32] i. Small initial infections die out asymp- 
totically as i -^- cxD. 

Proof: Since the infection strength satisfies ( |37] ). the 
matrix {f3^„A — SI) is Hurwitz. According to the property 
(|38T l of Gi(p,r), the system 

p = (/3,,A-«)p + Gi(p,0), 

which is system dSOl l with r = 0, is exponentially stable. In 
addition, according to Lemma [T] fi^-^j^,aijhj{Y) -^ oo 
as i ^' cx). Therefore, the term ri ^ a,ijPj in ( |26] l can be 
considered as a decaying disturbance for ( [30l l. Therefore, pi — ^ 
asymptotically as f — > oo. ■ 

Remark 6: From Theorem |2] the first epidemic threshold is 



(45) 



D. Epidemic Persistence in the Steady State 

The steady state is studied by letting the time derivatives pi 
and qi equal to zero, i.e., 

+ /3,gr 5]a.,pf -Jpf, (50) 

= K{i-pr - qn E «^^pf - ^-9^ E ^^^^f- ^^d 

From (ISTl i. it is inferred that 
Equivalently, according to (|52] i, the following is true 

ayPf. (53) 



cE«^.^r = ^E' 



Now, substitute for qf' ^ o.ijPT terms in ( |50] l using (|53T l 
to get 

/3o-^(i-K^)E«*^-pr 



p{Ay 



Pa 



1 



which is equal to the epidemic threshold in the classic SIS 
epidemic network. If the infection rate fi^ is such that 



5 p{Ay 



(46) 



the ratio -^ can be larger or smaller than -At, depending on 
the value of /3q. Therefore, if ( |46] l holds. Theorem [3] suggests 
that there exists another epidemic threshold r^. Using the 
definition of (3 in (l32T i. the condition ( |37] | in Theorem [3] 
can be expressed as 



P 



eq_ 
S 



B ^ 



Pa 



1 



^1 + ^ 

which is equivalent to 



< 



1 



5 l+P^- p{A) 



(47) 



\ 1 + % i + if/ 

(54) 

Theorem 4: Consider the SAIS epidemic model (fT3T l and 
(fT4l l. The steady state values of the infection probabilities of 
each individual in the SAIS model is similar to those of the N- 
intertwined SIS epidemic model ^ with a reduced infection 
rate /3^^. 

Proof: Based on the definition of /3g in (|32] |. the equation 
is simplified to 

/?e,(i-pr)E°^^pf-^pr = o, 

which can be expressed as 
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The second epidemic threshold r^ can now be obtained from 
inequality ( |48] | as 
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Notice that, according to ( |46] |. t\> t\. 
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Comparing (l55T l with dHJ from the Proposition |2] it is ob- 
served that the steady state values of the infection probabilities 
in an SAIS epidemic network is similar to those of a SIS 
epidemic network with reduced infection rate (3 . ■ 

Remark 7: The expression ( |32| | for fi can be rewritten as 

1 + lf 
The above expression suggests that (3 is always less than /3q 
since (3^ < (3^. It is insightful to look at the extreme cases 
for the values of (3 . Particularly, when the alerting rate k is 
very small, (3 -^ (3q, indicating that alertness plays a trivial 
role in the epidemic spread dynamics. When the alerting rate 
is very large, the reduced infection rate (3 -^ (3^. Another 



case, which is more important from the epidemiology point of 
view, is that if (3^ is very small, the epidemic spread can be 
completely controlled. 

V. Simulation Results 

Three examples are provided in this section. In all of the 
simulations, the curing rate is fixed at 5 = 1 so that the 
dimensionless time t — St is the same as the simulation time. 

Example 1: Consider a contact graph as represented in 
Fig. |2] For this network, the spectral radius is found to be 
p{A) = 3.1385. The alerting rate is arbitrarily selected as 
K = 0.1. The infection rate of an alert individual /3^ is chosen 
/3^ = 0.1. For the simulation purpose, nodes 1, 5, and 10 
are initially in the infected state. Other nodes are initialized 
in the susceptible state. In each simulation, the total infection 
fraction p{t) = j^^i=iPi[i) is computed. In Fig. [3] three 
trajectories are plotted. The trajectory (a) corresponds to the 
N-intertwined SIS model, with /3q = 2. Trajectory (b) is the 
solution of the SAIS model (fTSI l and (O developed in Section 
Hm Trajectory (c) is the solution of the SIS model but with 
the reduced infection rate (3 defined in (l32T i. As is expected 
from Theorem [T] the infected fraction in the SAIS model is 
always less than that of the SIS model. In addition, as proved 
in Theorem m the steady state infection fraction in the SAIS 
in equal to that of the SIS model with the reduced infection 
rate /3g^. 




Fig. 2. The contact graph in Example [T]and Example |2] 

Example 2: In Fig. [3] it can be observed that in the SAIS 
model the infection spreads similar to the SIS model at the 
first stage. Then, the size of the epidemics is reduced due 
to increased alertness in the network. In this example, for 
the same network in the previous example, the steady state 
value of the infected fraction and the maximum value of the 
infected fraction are presented as a function of the infection 
strength t = Pq/5. The simulation parameters are chosen as 
K = l, 13^ = 0.1. Note that 13^/5 = 0.1 < l/p^A) = 0.3186. 
Therefore, as discussed in Remark |6l there exists two distinct 
thresholds t], and r^ presented in ( |45] ) and ( |49] |, respectively. 
Simulation results for this example are shown in Fig. 2] 

Example 3: As is observed in Fig. HI the steady state values 
of the infected fraction p is zero before the second epidemic 
threshold r^. In addition, the maximum of the infected fraction 
is equal to the initial infected fraction before rj. The reason 
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Fig. 3. The infected population fraction in Example [T] (a) SIS model, (b) 
SAIS model, (c) SIS model with reduced infection rate fi^ 




Infection Strenglit r 



Fig. 4. The maximum infected fraction (blue hne) and the steady state value 
for the infected fraction (red line) in Example |2] 



for this observation is that before the first threshold t^, the 
epidemics dies out exponentially; as stated in Theorem |2] 
Between the two thresholds, maxtp(i) is greater than p(0) 
but steady state value p"* = 0. In other words, in this region 
the epidemic spreads at the first stage but then is completely 
controlled as a result of increased alertness. After the second 
threshold, p*^" < inaxtp{t), i.e., alertness reduced the size of 
the epidemic. 

Example 4: Consider an epidemic network where the con- 
tact graph is an Erdos-Reyni random graph with N — 320 
nodes and connection probability p = 0.2. The initial infected 
population is %2 of the whole population. The simulation 
parameters are /^g — 0.03, k — 0.05. Three trajectories are 
presented in Fig.|5] The trajectory (a) is for the SIS model, i.e., 
no alertness exists. Trajectory (b) is for /3^ = 0.02. In this case, 
the epidemic size is reduced in the steady state. Trajectory 
(c) corresponds to /3^ = 0.01, for which the epidemic dies 
out asymptotically. For the sake of evaluating the model 
development in Section |III] a Monte-Carlo simulation is also 
provided for each trajectory and shown in the figure in blue. 
As can be seen, there is a reasonable agreement between the 
proposed model (fTsT l and (fl4] i and the Markov process (|9|l. 
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Fig. 5. The infected population fraction in Example |4] (a) SIS model, 
(b) SAIS model with /3„ = 0.02. (c) SAIS model with (5^ = 0.01. The 
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VII. Conclusion 

In this paper, we add a new compartment to the classic SIS 
model to account for human response to epidemic spread. Each 
individual can be infected, susceptible, or alert. Susceptible 
individuals can become alert with an alerting rate if infected 
individuals exist in their neighborhood. An individual in the 
alert state is less probable to become infected than an individ- 
ual in the susceptible state; due to a newly adopted cautious 
behavior The problem is formulated as a continuous time 
Markov process on a general static graph and then modeled 
into a set of ordinary differential equations using mean field 
approximation method and the corresponding Kolmogorov 
forward equations. The model is then studied using results 
from algebraic graph theory and center manifold theorem. 
We analytically show that our model exhibits two distinct 
thresholds in the dynamics of epidemic spread. Below the first 
threshold, infection dies out exponentially. Beyond the second 
threshold, infection persists in the steady state. Between the 
two thresholds, the infection spreads at the first stage but then 
dies out asymptotically as the result of increased alertness in 
the network. Finally, simulations are provided to support our 
findings. Our results suggest that alertness can be considered 
as a strategy of controlling the epidemics which propose 
multiple potential areas of applications, from infectious dis- 
eases mitigations to malware impact reduction. Generalizing 
the current results to time-varying weighted topologies is a 
promising extension. 
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